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Relevance of intra- and inter-subband scattering on the absorption in 
heterostructures 
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We analyze the absorption lineshape for inter-subband transitions in disordered quasi 2D heterostructures by 
an exact calculation. The intra-subband scatterings control the central peak while the tails of the absorption 
line arc dominated by the inter-subband scattering terms. Our numerical study quantitatively assesses the 
magnitude of the free carrier absorption. The accuracy of different models currently used for gain/absorption 
is discussed. 
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The demand to produce reliable terahertz and infrared 
devices has triggered extensive research on intersubband 
transitions in semiconductor heterostructures^. A high- 
light is the realization and continuous improvement of 
Quantum Cascade Lasers (QCL)^"^ that emit in this 
frequency range. The search for improved performances 
of QCLs includes a better design of the layer sequence 
to enhance the population inversion between the two 
subbands involved in the lasing transition, the decrease 
of the non radiative paths, but also the control of 
the photon re-absorption by the free carriers that are 
present in the structure^^^^. The free carrier absorption 
(FCA), well known from bulk structures, is transferred 
to a variety of intra-subband and inter-subband oblique 
(in 2D k space) transitions in heterostructures^^'^"'. 
Tailoring these (often parasitic) processes requires a 
thorough understanding of the nature of the optical 
transitions in imperfect heterostructures. 

For light polarized along the growth (z) axis of an ideal 
semiconductor heterostructure, the energy dependence of 
the absorption coefScient between subbands 1 and 2 is a 
Dirac delta function line centered at Hluq = E2 — Ei^^. 
This "atomic-like" profile stems from the fact that the 
electromagnetic wave coupling requires a non-vanishing 
z-dipole (1 1^x12), but does not affect the in-plane motion 
of the electrons. These motions are identical in the two 
subbands (i.e. plane waves, owing to the in-plane trans- 
lation invariance) and their corresponding energies form 
identical dispersions (if we assume the same parabolic 
mass TO*). Static disorder affects the in-plane motion 
by introducing intra-subband (V^^f"^^) as well as inter- 
subband (V'J^f'''') couplings. Bound states usually ap- 
pear below the 2D subband edges. As a consequence, 
the strict selection rules leading to the Dirac lineshape 
are no longer expected to apply to a disordered structure 
which is frequently taking into account by some effec- 
tive broadening. Here we treat this feature from a fully 
microscopic point of view by considering the exact eigen- 
states of the disordered quantum well potential. Our 



calculations allow us to gain a complete insight into the 
underlying physical processes leading to the absorption 
in heterostructures, not only near resonance, but also far 
above or far below the resonance (theses regions have not 
deserved enough attention in the literature, even though 
of great technological interest). We will show that the re- 
gions far from resonance can nicely be understood by re- 
curring to defect-induced scattering processes, and high- 
light the respective roles of VJ"*''^ and VJ^j Finally, the 
complete calculation will allow us to determine the accu- 
racy of different models currently employed to evaluate 
the gain/absorption of such structures. 

In the presence of disorder, the heterostructure cigcn- 
problem can be written as 'i!^{p,z) = F^.n{p)Xn{z), 
with: 



2to* 



n' 

= (£,-£;„) F,,„(p) (1) 



where p={x,y), F^^n{p) is the in-plane envelope 
function of the n*'* subband and V^"'" {p) = 
{Xn{z)\Vaef{p, z)\xn'{z))- Xn{z) IS the heterostructurc lo- 
calized part and En the corresponding bound state. In 
absence of defects v {n,k), Fy[p) — ^ exp(zA: • p)/y/S 
and Ely — > En + h^k"^ /{2m*). In the following we shaU 
present results obtained after a numerical diagonaliza- 
tion of Eq. (1) within a truncated basis with the two 
lowest subband states (n=l,2). In practice, we expand 
Fv,n=i.2{p) in a plane wave basis that fulfill periodic 
conditions in a 200x200 nm^ box. The matrix cle- 
ment of the dipolc coupling to light are proportional 
to {pz){u,i_i) = {^y{r)\pz\^ ii{r))- For defect-free struc- 
tures one has the selection rules: {Pz)(^n k) («' t?) ~ ^ 
Sn,n')Sj: g{xn\Pz\Xn'), reflecting the in-plane translation 
invariance and the non-vanishing z-dipole. In the follow- 
ing, we will be interested in the effect of disorder on these 
matrix elements. 

In QCLs interface defects constitute an important clas- 
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tic scattering mechanisms, which we will consider in de- 
tail here. Most features can be directly transferred to 
other mechanisms such as impurity and alloy scatter- 
ing. Interface defects arise from the non ideality of the 
well/barrier interface at z=zq between two consecutive 
layers of the heterostructure: they can be a protrusion of 
the barrier material in the well (repulsive defects) or vice- 
versa (attractive defects). For nearly ideal interfaces, the 
interface defects are only one monolayer thick. The scat- 
tering potential created by one disordered interface has 
the form 

^def(r) = f{z)udel{p} (2) 

showing a separation on the dependence on the z and p 
variables. Udef{p) is the superposition of localized func- 
tions centered at the randomly placed interface defects on 
the (x, y) plane. For a barrier/well interface at z = zq and 
for a repulsive defect f^°^{z) = VbQ{z — zo)Q{zq + h~ z) 
while for an attractive defect f^*'^{z) = —Vb'd{zo — 
z)'d{z — zq + h), where Vf, is the barrier height, h is the 
thickness of one monolayer, Q{z) is the Heaviside func- 
tion. 

In Fig. 1 we show the matrix |(^',y(r)|p2|5'p(r))p for 
our calculated eigenstates of the disordered heterostruc- 
ture. Here we consider a 9/2/3 nm GaAs/Gao.75Alo.25As 
double quantum well (DQW) structure with Gaussian in- 
terface defects placed on the two inner interfaces, with a 
characteristic in-plane extension of 3.6 nm and a frac- 
tional coverage of the surface N^dT^a'^ / S = 30%. The 
figure clearly displays two blurred straight regions around 
\£fj. — Ei/I ~ ftwo = 73.8 meV, corresponding to the sub- 
band spacing in this sample. If there where no disorder, 
there would be no blurring since (see above) a single final 
state would match any given initial state. The fact that, 
the matrix element is almost zero if the energy difference 
between the true states differs strongly from the intcr- 
subband spacing corresponds well with the conventional 
broadening picture. 

In this context it has to be noted, that the blurring 
is strongly reduced, if the wave functions Xniz) for the 
subbands n = 1 and n — 2 arc of the same magnitude 
at each interface. Neglecting V|]"f°'' this provides iden- 
tical (or very similar) in-plane wave functions for both 
subbands, which provide strong selection rules for the z- 
matrix elements. Similar selection rules appear in strong 
magnetic fields^^. In our case, the wave functions differ 
essentially at the interfaces, as the excited state pene- 
trates stronger into the barrier, so that this effect is only 
of very minor relevance. 

We show in Fig. 2 (red full line) the absorption coef- 
ficient of the DQW structure based on the matrix ele- 
ments of Fig. 1 (in all absorption calculations, the delta 
of energy conservation was replaced by a Gaussian with 
a FWHM=1.88 meV). In this absorption spectrum, the 
contributions of V"J"f'''' and V'drf'''^ ^'^^ taken at all or- 
ders and are fully admixed. In order to highlight the 
physics underlying the assisted light absorption and dis- 
entangle the various contributions to its spectrum, we 
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FIG. 1. Color plot of the decimal logarithm of the normal- 
ized squared modulus of the optical matrix elements for the 
ensemble of transitions e^, — in a disordered 9/2/3 nm 
GaAs/Gao.75Alo.25As double quantum well. The disorder is 
due to interface defects randomly distributed on the two inner 
interfaces of the structure (see text). 



discuss in the following the results obtained within dif- 
ferent approximations and/or models. We also show in 
Fig. 2 the calculated absorptions when only V^™* (green 
dotted line) or only Vj^f'^'^ (orange dashed-dotted line) 
are retained in the scattering matrix. We see that the 
low energy (the high energy) FCA is dominated by the 
Kiof°'^ i^dci'^^) couphngs. Indeed, the initial and final 
states involved in the non-resonant absorption displays a 
dominant Ei or E2 subband character. Thus, since the 
dipole coupling only triggers inter-subband events and 
starting from an initial state in Ei, an additional inter- 
subband scattering is needed for the low energy intra- 
Ei oblique transition. Conversely, the high-energy ab- 
sorption tail involves essentially Ei-to-E2 processes and 
thus rely mostly on V^c^^- is worth pointing out that 
the near-resonant absorption (see Fig. 2b) is accurately 
rendered by the only-intra approximation, whereas the 
only-inter one fails in predicting both the intensity and 
the absorption profile. Finally, we also show in Fig. 2 the 
perturbative estimates of the FCA obtained by expand- 
ing the electron wavefunction to the first order in Vdcf^^ 
(dashed blue line). We see that the perturbative estimate 
provides an excellent rendering of the shape and strength 
of the actual absorption when it is justified, i. e. not too 
close to the resonance. 

An important issue is to understand how the results of 
the exact diagonalization compare to correlation function 
approaches, essentially based on Green's function meth- 
ods. To this end, we show in Fig. 3 the absorption coef- 
ficient for the El — E2 transition of the DQW structure 
described above, calculated with two different correla- 
tion function (CRF) approaches, and we compare it with 
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FIG. 2. Absorption spectrum for the E1-E2 transition calcu- 
lated by; fully numerical diagonalization (red full line); tak- 
ing into account either only Vj"f (green dotted line) or only 
Vj"f '^'^ (orange dashed-dotted line) ; and by expanding the elec- 
tron wavefunction to the first order in both V^^f^^ and Vjof'^"^^* 
(blue dashed line). Inset: Zoom of the absorption spectrum 
in (a) around tiuo. T = 100 K. 



the one calculated by exact diagonalization. In CRFi 
the Keldysh Green's function formalism was applied us- 
ing the implementation of Refs. 17 and 18. In order to 
make a quantitatively meaningful comparison among the 
various calculations, we numerically extracted the cor- 
relation length and the average defect depth from the 
randomly generated interface defects used in the exact 
calculation without the use of any fitting parameters. We 
see that CRFi gives a good description of the low uj be- 
havior, as the inclusion of the nondiagonal self-energies^^ 
fully covers the inter-subband terms addressed above. At 
large energy, CRFi overestimates the absorption; a fea- 
ture that is related to the use of fc-independcnt scattering 
matrix elements in order to simplify this numerics. This 
assumption of an effective delta- potential scattering po- 
tential overestimates the scattering for large wave- vector 
transfer, which is crucial for the tails. For the peak, 
however, this approximation turns out to provide a full 
quantitative description, without any adjustable param- 
eters. In CRF2, we use Unuma et al.^^ formalism for the 
intersubband absorption coefficient that follows Ando's 
approach^*^. To draw CRF2 we converted Unuma et a/.'s 
real part of the conductivity into an absorption coefficient 
following^^ and we used the interface disorder potential 
described above. It is clear that CRF2 poorly describes 
the far from resonance FCA. This is most probably due to 
the fact that, contrarily to CRFi, the off-diagonal com- 
ponents of the Green's function are neglected, whereas 
the diagonal ones are evaluated within the self-consistent 
Born approximation. 

In conclusion, we have computed numerically the line- 
shape of inter-subband transitions in hcterostructures 
with interface disorder. We have found that only few fi- 
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FIG. 3. Absorption spectrum for the E1-E2 transition calcu- 
lated by three different models: exact diagonalization (red full 
line), Keldysh Green's function formalism, CRFi^^'^* (blue 
dashed line), Unumas model, CRF2^^ (green dotted line). 
T = 100 K. 



nal states are optically connected to a given initial state 
despite disorder, which reflects the common picture of 
a broadened transition. While intrasubband scattering 
dominates at the resonance peak and at higher photon 
energies, intersubband scattering dominates the low en- 
ergy tail of the absorption spectrum. We have compared 
the outputs of several models to compute absorption to 
the fully numerical results, in particular their low w be- 
haviors. The perturbative estimate for FCA^^ works 
very well for the absorption tails. We find also that 
Unuma et aZ.'s model"'^^ gives a poor description of the 
non-rcsonant FCA tails, while the Keldysh Green's func- 
tion formalism^^ with full nondiagonal self-energies fits 
the numerical results nicely. A clear picture of the inter- 
subband and intra-subband transitions now emerges from 
this study and emphasizes the need for a proper account 
of the perturbation of the carrier wavcfunctions by de- 
fects and not only of their energies. 
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